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Abstract 

Structure in the Universe grew through gravitational instability from very smooth initial condi- 
tions. Energy conservation requires that the growing negative potential energy of these structures 
is balanced by an increase in kinetic energy A fraction of this is converted into heat in the col- 
lisional gas of the intergalactic medium. Using a toy model of gravitational heating we attempt 
to link the growth of structure in the Universe and the average temperature of this gas. We find 
that the gas is rapidly heated from collapsing structures at around z ~ 10, reaching a temperature 
> 10 6 K today, depending on some assumptions of our simplified model. Before that there was 
a cold era from z ~ 100 to ~ 10 in which the matter temperature is below that of the Cosmic 
Microwave Background. 
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Introduction. — It is well known that the temperature of the Cosmic Microwave Background 
is T 7 o = 2.725 K, while the atomic matter today is typically orders of magnitude hotter than 
this. As the Universe expands the photon temperature drops so that T(z) = T 7 o(l + z), 
and Compton scattering couples the matter to the radiation until about z — 100, when it is 
allowed to cool adiabatically with T oc (1 + z) 2 . In a perfectly smooth Universe with the 
same background cosmology as we observe today, the present temperature of the gas would 
be only about 20 mK [l|. 

However, the Intergalactic Medium (IGM) is composed of gas at a variety of temperatures 
(e.g. 0]), ranging from ~ 10 4 K in dense clouds, up to > 10 6 K in the rarest inter-cloud gas. 
The reason that the gas is so hot is patently because the Universe contains structure. Our 
lumpy Universe generates thermal energy in a number of ways, but the simplest source is 
purely gravitational. The importance of gravitational heating for understanding the struc- 
ture and evolution of galaxy clusters has been discussed by a number of authors (e.g. see [!, 0] 
for recent examples). However, an explicit connection between gravity and heating of the 
IGM seems to have escaped notice in the literature. The gas is effectively heated from the 
growth of structure through gravitational instability, with the thermal energy coming from 
the increasingly negative energy of the growing potential wells. 

The idea of gravitational shock heating of the IGM as a direct result of structure formation 
goes back at least to the 'pancake' model of the 1970s This work was later extended 
using the 'Zel'dovich approximation' to follow the shock-heating of gas outside collapsed 
objects (see e.g. [6|), or in a related approach to use an extension of the Press-Schechter 
formalism to estimate the fraction of shocked gas Q. Such numerical calculations allow for 
an investigation of the contributions of collapsed and shocked gas to temperature evolution 
of the IGM (see also 0, B E3)- 

The details will of course be quite complicated, not least because starbursts and quasars 
provide photon and mechanical heating to the IGM in very non-linear and inhomogeneous 
processes. However, we will here focus only on gravity and aim to estimate the temperature 
of the IGM from the potential energy of the Universe, using a toy model to highlight the 
basic connection between gravitational and thermal energies. For our numerical work, we 
use the baryon density £\ = 0.045, cold dark matter density Q c = 0.0255, Hubble parameter 
H = 70kmsec _1 Mpc -1 , and the spectral index of initial fluctuations n s = 1. To normalize 
the power spectrum we fix the matter variance in 8 h~ x Mpc spheres to cig = 0.9. These 



values are consistent with the current best fit cosmological parameters [11 



Potential energy estimate from power spectrum. — The volume averaged gravitational poten 



tial energy (GPE) per unit mass as a function of redshift z is [12 



^) = |((i + ^zMv)), (i) 

where <fi is the Newtonian potential, defined by the line element ds 2 = 

a 2 [— (1 + 2(f)) dr + (1 — 2(p)dx 2 ], and r is conformal time. Using the Poisson equation 
V 2 </> = 4irGa 2 p5 and the definition of the correlation function 

£(r, z) = (5(x, z)5(x + r, z)) = JL J d *kP(k, z)e^ r , (2) 
where P(k,z) = \6(k, z)\ 2 , the Newtonian potential energy is 

W(z) = - f ^( k,z) , (3) 
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FIG. 1: Normalized gravitational potential as a function of scale at z = (solid), 1 (dotted), 10 
(short-dashed) and 100 (long-dashed). The corresponding curves with additional small-scale power 



are non-linear corrections from the HAL0FIT code 141 ] . 



and with the dimensionless gravitational power 



Al,(k,z) = 3H ^f + z) P(k,z)k. (4) 



This can also be obtained from the definition in Peebles [l3[], W = — ^Gp m a 2 J d 3 r £(r, z)/r = 
— 2nGp m J 2 , where J 2 = J dkP(k, z)/(2ir 2 ). In Fig. [T] we show (1 + z)A 2 N {k,z) for our 
cosmological model at z = 0, 1, 10 and 100. The peak of the gravitational energy contribution 
occurs at a comoving scale of ~ 10 Mpc. 

In the matter dominated era (z > 1) the quantity (1 + z)A 2 N (k, z) remains roughly con- 
stant, since the growth factor, given by D(z) = S m (z)/5 m (z = 0) , scales as (1 + z)^ 1 . 
At lower z, in the dark energy dominated era, the growth of structure is significantly sup- 
pressed due to the increasing expansion rate - this slow down has recently been detected in 



observations of galaxy clusters [15 



In order to estimate the IGM temperature T gas (z) one can suppose that the GPE is 
equal to the average kinetic density energy of baryonic matter. Therefore, one could set 
p m W = pbkT/mn, where k is Boltzmann's constant, mn the mass of hydrogen (neglecting 
helium for simplicity) and p m , pb are the energy densities of matter and baryons respectively. 
We show the results of this calculation in Fig. [2J along with T(z) from the recombination 



code RECFAST [161. 1171]. which assumes the matter distribution is perfectly smooth and cools 
adiabatically. One finds an average temperature of ~ 10 8 K at z = 0, and the 'gravitational 
temperature' exceeds the RECFAST value for z < 1000. 

Clearly, this linear calculation over-predicts the temperature. The reason for this is that 
gravitational shock-heating is associated with the collapse of non-linear objects and shell- 
crossing. The linear GPE is largest on scales of ~ 10 Mpc, which are only just going 
non- linear at z = 0. Linear scales are associated with smooth bulk flows, and hence there 
is no mechanism for the gas to be heated. To be more realistic, we should only include 
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FIG. 2: Evolution of T gas (z) for linear theory (dotted curve), non-linear theory with a cut-off scale 
(solid) and halo number density (short-dashed) . We also show the baryonic temperature evolution 
from the RECFAST code lall3] (long-dashed) and the CMB temperature (dot-dashed). 



the total GPE coming from non-linear scales at a given redshift. In this regime, setting 
p m W = pbkT/mii implicitly assumes that the damped outgoing shock waves from collapsed 
objects efficiently (and instantaneously) heat the IGM. 

Since the linear power spectrum becomes becomes inaccurate on small-scales, we include 
non-linear corrections to P(k, z) from the the HALOFIT code 
shown in Fig. (TJ and result in an increase of small-scale power. 
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These corrections are 
In order to estimate the 
cut-off scale at which perturbations are going non-linear, we compute the value i?NL at which 
the RMS mass variance 



4( 



j ^A\k,z)W 2 (kR) 



(5) 



is equal to unity. Here W(x) is the window function associated with a spherical top-hat, 
and A 2 (/c, z) = k 3 P(k, z)/(2n 2 ) is the dimensionless power spectrum. We then perform the 
integral only above &nl = 27r/i? NL . The results of this computation are shown in Fig. [2]- 
one finds a smaller temperature at z = 0, and a much faster decrease with redshift. This is 
due to /cnl shifting to smaller-scales for increasing z, so less GPE contributes to the heating. 
Halo number density estimate. — We can also estimate the gravitational energy in a differ- 
ent way by just considering viralized objects. Virialization occurs at higher over-densities 
than those discussed previously - in a flat Einstein-de-Sitter model linear theory predicts 
a spherical top-hat will break away from the expansion at a mass fluctuation of 5 = 1.063, 
collapse at 5 = 1.686, but virialize at 5 = 1.59 (see e.g. [3]). We use the Press-Schechter 
mass function 19[ to estimate the number density of viralized objects 



dn 
dM 



z,M) 



2 p m (t 



M 2 



dhiay[{z = 0) 
d\nR 



(6) 



where v = 5 C / [D(z)ayi{z = 0)], the critical threshold 5 C = 1.686, the mass enclosing a sphere 
of radius R is M = 4:iip m {t )R 3 /3 and the growth function is normalized to unity today. 
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The growth function can be computed from 



5fi m ,o H(a) f a da' 



(7) 



2 H J (a'H(a')/H ) 
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where the scale-factor a = 1/(1 + z) and H(a) is the Hubble rate. 

In order to compute the gravitational energy we can use the virial theorem, W = —2K, 
where the kinetic energy per unit volume is K — k{T V i r )M toi . The mass averaged virial 
temperature is defined as 



with M tot = fp°° M-^jjdM. We use the virial temperature normalization for a conventional 
cosmology in H, with T vir (K) = 5 x 10" 23 (M/kg) 2 / 3 (l + z). 

We again assume that outgoing shock waves effectively share out the GPE among all 
the particles in the Universe, not just those in virialized structures. Hence one can equate 
k(T vir )M tot = phkT to obtain the evolution of T(z). This is shown in Fig. [2]- the redshift 
evolution has a similar profile to the non-linear power spectrum estimate (with a cut-off 
scale -Rnl); although the temperature is approximately an order of magnitude lower. One 
can see that effectively these two calculations are very similar, with mainly just a different 
choice of cut-off scale. 

Gas cooling. — So far, we have assumed that the IGM is heated without taking into account 
cooling processes in the gas. We can expect this to be a reasonable estimate as long as the 
cooling time-scale is longer than the Hubble time H~ x . Assuming the cooling is dominated 
by thermal bremsstrahlung we show the ratio of cooling to Hubble time in Fig. [3J For our 
model of IGM heating from only collapsed objects, we find the cooling time becomes less 
than the Hubble time around z ~ 10. At this point the baryonic matter in the IGM will 
begin to heat from its adiabatically cooled value of ~ 3 K. The CMB temperature is ~ 30 K 
at z ~ 10, so it appears there is an epoch from z ~ 10-100 when the matter is actually colder. 
Since the IGM is outside collapsed regions and thus taking part in the Hubble expansion, it 
will also continue to cool adiabatically - including this additional effect in our computations 
leads to a reduction in temperature by a factor of 2-3 at z = 0. 

Conclusions. — Since the Universe is lumpy, it is necessarily hot. A simplistic picture is that 
the temperature of the IGM comes from energy balance with a fraction of the gravitational 
energy that is building up through gravitational instability. Virliaization of structures re- 
quire that about half of the potential energy is lost to the 'environment', here meaning that 
the dark and baryonic matter acquire kinetic energy. When shell-crossing occurs kinetic 
energy is converted into thermal energy in the collisional material, the IGM gas. We have 
shown that this picture leads to an IGM which is significantly colder than the CMB until 
gravitational heating takes over at z~ 10. 

Things became rapidly more complicated just after this era. An accurate calculation of 
the heating process would require solving for the inhomogeneous growth of structure, in- 
cluding hydrodynamic effects, as well as cooling processes. In addition, reionization of the 
Universe appears also to happen at z ~ 10 and so photon sources, electromagnetic inter- 
actions and radiative transfer need to be considered in order to fully understand the IGM 
today. Many astrophysics theorists are working hard on just these problems. 
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FIG. 3: Ratio of cooling time to Hubble time. Lines are labelled the same as in Fig. [2j 

discussions on this topic, in particular Alan Duffy, Kris Sigurdson, Martin White and Jim 
Zibin. 
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